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Abstract 

We calculate the dominant O [G'pMfj) two-loop electroweak corrections to 
the fermionic decay widths of a heavy Higgs boson in the Standard Model. 
Use of the Goldstone-boson equivalence theorem reduces the problem to one 
involving only the physical Higgs boson H and the Goldstone bosons 
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and z of the unbroken theory. The two-loop corrections are opposite in sign 

to the one-loop electroweak corrections, exceed the one-loop corrections in 

magnitude for Mji > 1114 GeV, and increase in relative magnitude as M^j 

for larger values of Mh- We conclude that the perturbation expansion in 

powers of GpHfj breaks down for Mh ~ 1100 GeV. We discuss briefly the 

QCD and the complete one-loop electroweak corrections to H —>■ bb, tt, and 

comment on the validity of the equivalence theorem. Finally we note how a 

very heavy Higgs boson could be described in a phenomenological manner. 
PACS number(s): 12.15.Lk, 11.15.Bt, 14.80.Bn 
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Typeset using REVT^ 



I. INTRODUCTION 



One of the great puzzles of contemporary elementary particle research is whether nature 
makes use of the Higgs mechanism to generate the observed particle masses. In the minimal 
standard model (SM) of electroweak interactions, the symmetry breaking is implemented 
using this mechanism with one weak-isospin doublet of complex scalar fields with weak 
hypercharge Y = 1. With the spontaneous breaking of the SU(2)i(8)U(l)y gauge symmetry, 
three of the four scalar degrees of freedom are absorbed to create the longitudinal polarization 
states of the intermediate bosons, and Z. At the same time, the quarks and charged 
leptons acquire masses through their Yukawa interaction with the scalar doublet. There 
remains at the end one neutral scalar boson with positive parity and charge conjugation, 
the physical Higgs boson H. 

Most of the properties of the scalar or Higgs sector of the SM are fixed experimentally, 
e.g., the vacuum expectation value, v = ~ 246 GeV, the coupling of the Higgs 

boson to the gauge bosons, gvvH = 

25/4(^1/2^2 ^i^g^g y ^ w^,Z, and the coupling of 
the Higgs to the fermions, gjfH = 2^/^G]/^m/. However, the mass Mh of the Higgs boson 
and its quartic self-coupling, A = GpMfj / \/2^ are unspecified. It is therefore of considerable 
interest to analyze processes which can give theoretical limits on Mh-, or test the effects of 
the quartic coupling phenomenologically. 

The range of possible Higgs masses is constrained from below both experimentally and 
theoretically. The non-detection of the Z-boson decay Z ffH at LEP 1 and SLC has 
ruled out a Higgs mass of less than 63.9 GeV at the 95% confidence level Depending on 
the mass of the top quark, the requirement that the vacuum be the true ground state could 
provide an even more stringent theoretical lower bound 0. Other theoretical arguments 
bound the Higgs mass from above. Nonperturbative lattice computations give an upper 
limit for Mh of about 710 GeV [Q. Unitarity arguments in intermediate-boson scattering at 
high energies and considerations concerning the range of validity of perturbation theory 
I^J^ establish an upper bound < (87r\/2/3G'i?)^/^ ~ 1 TeV in a weakly interacting 
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theory. The unitarity bound on Mh is lowered significantly when the approach of [^,0 is 
extended to higher orders; see and references therein. However, the improved bound 

depends on the energy scale up to which the SM is assumed to remain valid. 

A violation of the unitarity bound on Mh is presumably a signal for the onset of strong 
interactions in the Higgs sector of the SM, a possibility which is of considerable interest in 
its own right It would therefore be desirable to sharpen the bound by removing the 



uncertainty associated with the mass scale at which it is applied, or to find a separate scale- 
independent bound. In fact, the work presented here on two-loop electroweak corrections 
to the fermionic decay modes of the Higgs boson, H —>■ ff, gives a scale-independent 
limit on in a weakly interacting theory. We find that the two-loop corrections to the 
fermionic decay rates exceed the one- loop corrections in magnitude for ^ 1114 GeV, 
and increase in relative magnitude proportionally to Mjj for larger Higgs-boson masses. We 
conclude, as reported previously [p!2|-[l^, that the perturbative expansion fails to converge 
satisfactorily, and that the theory becomes effectively strongly interacting in the Higgs sector 
for Mh > 1100 GeV. 

The result noted above is a consequence of our calculation of the dominant two-loop 
electroweak corrections to the fermionic decay rates of the Higgs boson, the subject of 
this paper. The fermionic decay modes are of considerable phenomenological interest. For 
example, the Higgs boson decays predominantly to bb pairs if Mh ^ 135 GeV. The search 



for a low- or intermediate-mass Higgs boson at future high-energy e~^e linear colliders [ITS 



the Fermilab Tevatron [0, or a possible 4-TeV upgrade thereof JT^ will rely largely on this 
mode by tagging the B mesons. Moreover, in this mass range, the branching fractions of all 
other decay channels depend sensitively on the H ^ bb decay width. It has been argued that 
the low-mass Higgs boson might also be detectable at future hadron supercolliders through 



the H r^r^ signal |T8[, while the decay H ^ tt will have an appreciable branching 
fraction for Mh > 2mt. Future high-energy e"'"e~ colliders will also be able to measure the 
Hti Yukawa coupling |T^ . 



The measurement of the Higgs-boson mass and couplings in future experiments will 



require an understanding of the radiative corrections to the fermionic decay rates of the 
Higgs boson. Much work has been done in this area, and a recent review is given in pO[ . 
There are important differences between the radiative corrections involved in Higgs physics, 
and those famihar in the gauge sector of the SM, e.g., in Z-boson decays. It is well known 
that the quantum effects induced by virtual Higgs bosons are screened in Z-boson physics 
||7|j2l||: they depend only logarithmically on Mh at one loop, and are quadratic in Mh-, but 
with minute coefficients, at two loops In contrast, the one-loop electroweak corrections 
to the partial decay widths |f7ip3H25|| and production cross sections of the Higgs boson 
are already dominated for Mh ^ My/ by terms quadratic in Mh- These terms give rise 
to moderate enhancements of the rates for Higgs masses of up to 1 TeV. However, it is 
premature to conclude that the two-loop electroweak corrections will also be perturbatively 
small in the YA^-Mh range, since these corrections have terms quartic in Mh- It is of both 
theoretical and phenomenological interest to check the importance of these potentially large 
corrections by explicit calculation. 

In this paper, we calculate the complete O {G']^Mjj) corrections to the fermionic decay 
rates of a Higgs boson with Mh 3> M\y. These corrections, which are the leading two- loop 
electroweak corrections for Mh ^ My/, are independent of the fermion flavor, and, as noted 
above, are larger than the one-loop corrections of O {GpM'jj) for Mh > 1114 GeV. We 
compare our results with other known one-loop corrections. To obtain the full electroweak 
two-loop corrections for specific fermion channels in the limit Mh ^ M^, one would have 
to calculate further, flavor-dependent corrections of mixed orders in the Higgs and Yukawa 
couplings, namely O (^G'j^'Mfjmj^ and O (^G'jpmj^ . These corrections, however, are not uni- 
versal. For example, different fermionic channels, such as H ^ t^t^, bb, ti, all have 
different dependence on m^. 
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II. CALCULATION OF THE O {Gl^Mfj) CORRECTIONS 

In this section, we sketch the calculation of the dominant, flavor- independent electroweak 
corrections to the decay rate H ff. The starting point of our analysis is the bare 
Lagrangian for the Higgs-fermion interaction, 

4o' = -^^/,o^o^/,o, (2.1) 

where the subscript "0" denotes bare quantities. Our aim is to obtain the leading, flavor- 
independent corrections to the Hff vertex in powers of GpMjj. The mass and wave- function 
renormalization constants for the fermions as well as the loop corrections to the vertex 
depend on the Yukawa couplings for the fermions, and are omitted consistently because 
of their subleading nature. (We will return to these corrections at the one-loop level in 
Sec. IIIB.) We may therefore replace the fermionic quantities mf^ and ipf^ in Eq. 
by m/ and ipf. The contributions to the renormalization constants for the bare Higgs 
field, Hq, and the bare vacuum expectation value, Vq, in powers of GpMfj are determined 
entirely by the symmetry-breaking sector of the SM. The contributions of fermion loops 
to these two quantities depend on the Yukawa couplings and are again omitted as they do 
not contribute to the 0{G\Mfj) corrections (see Sec. IIIB). We may therefore calculate the 
desired corrections to the decay H ^ ff vertex with all Yukawa couplings set to zero. 
The calculation can be simplified greatly in the limit of interest, ^ Myy, through the 



use of the Goldstone-boson equivalence theorem ||2^. This theorem states that the leading 
electroweak contribution to a graph in powers of GpM'jj can be calculated by replacing the 
gauge bosons W"^, Z by the would-be Goldstone bosons w^^ z of the symmetry-breaking 
sector of the theory. Because Mw/Mh oc qv/Mh, we can simplify our calculation consis- 
tently in the limit of a heavy Higgs boson by neglecting the gauge couplings g, g' and taking 
the Goldstone bosons to be massless. Adopting the conventions of p8[, we can write the 
relevant Lagrangian for the symmetry-breaking sector of the SM in terms of bare quantities 
as follows: 
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(2.2) 



A, 



+ Hl)^ - Xovo (wl + Hi) Ho , 



where the real scalar triplet, w = {wi,W2, w^), is related to the Goldstone bosons, and z, 
by = {wi^iw2) I \pl and z = W3, respectively. The tadpole counterterm, which cancels all 



tadpole contributions of jCq^^ order by order, has been omitted in writing Eq. 



therefore 



all graphs which include tadpole contributions need to be dropped in calculations p8 



Note that a full equivalence-theorem calculation would also require the complete Yukawa 
Lagrangian for the interactions of fermions with the massless Goldstone bosons and the 
Higgs boson [^]. Because we are not interested in corrections due to the Yukawa couplings 
we do not give the complete Yukawa Lagrangian here, except for the piece given in Eq. (2.1). 

The on-mass-shell renormalization is carried out in such a way that the physical mass of 
the Higgs boson, defined in terms of the position of the pole in the Higgs propagator, is M^- 
The three Goldstone bosons remain massless to all orders in the perturbation expansion in 
Ao, and satisfy a residual S0(3) symmetry |]30|. Requiring that the residues of the physical 
on-shell propagators be unity fixes the wave-function renormalization constants defined by 
the relations Wq = Z^J'^w^, zq = Z^^z, and Hq = Z]I^H, where w^, z, and H are the 



physical fields. The result is p8 



— = 1 - —no (p2^ 

1 



(2.3) 



_d_ 

dp'^ 



p2=0 



p2=0 



where n° (p^), n°(p^), and U'jj^p'^) are the self-energy functions for the bare fields calculated 
from the Lagrangian in Eq. (^^) . Because of the S0(3) symmetry, H° = H° and Z^ = Z^. 
Explicit expressions for the H°'s correct to two loops may be found in Eqs. (11) and (12) of 
112 811 . Furthermore, we have 1281 



Mlo = Ml- ReH°, (iVfi) , 



(2.4) 



Ml, = -ReUliO) = -UliO) , 

A / Ren", (M^)-^°(o) ^^ 

Using these results, we can write the Lagrangians above entirely in terms of renormalized, 
physical quantities. The physical vacuum expectation value is fixed in terms of the Fermi 
constant Gp by the familiar relation v = , while the physical quartic coupling 

is given by A = GrMfj/ \/2. The renormalized symmetry-breaking Lagrangian is given in 
Eq. (7) of ||28|, while the renormalized form of the Higgs-fermion Lagrangian is given for our 



purposes by 



<^r--%'-^i>^^!- (2-5) 



The radiatively corrected fermionic decay rate of the Higgs boson is consequently given by 

t{h^ ff) = ^Vb{h^ ff) , (2.6) 



where 31 



is the Born result. Here A^c = 1 (3) for lepton (quark) fiavors. 

The wave-function renormalization constants Zu and were calculated to two loops, 
O {G'jpMfj), in [^] using dimensional regularization. The two-loop diagrams that contribute 
to Zh and Z^j at that order through the derivatives of the self-energy functions in Eq. ( p.3|) 
are shown in Fig. [l|; no single diagram gives an exceptionally large contribution. The results 
of the calculation can be written in the form 

^ = 1 + Ar (a. + 0(e)) + X^e^ (- + 6. + 0(e)) + 0(A=^) {a = w,H), (2.8) 



where A = (A/IGtt^), ^ = Airfi"^ /Mfj, e = (4 — D)/2, D is the dimensionality of space-time, 
and /i is the arbitrary scale parameter introduced in the interaction to keep A dimensionless 
for e 7^ 0. The coefficients in the expansions above are: 
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= 1 , (2.9) 

bn, = ^ + 2C(2) - 67 - SnVs + UCh (^0 Vs 

^ 6.098, 
aH = -12 + 27rV3 

^-1.12, 

fen = ^ - 96C(2) + 90C(3) - 67 - 487ra2 (^) + IWnVs - 216012 (^) Vs - I62K5 



41.12 



The constant = 0.92363... was evaluated numerically from Eq. (A86) of pSf. The 
Riemann ( function takes the values ({2) = 7r^/6 and ({3) = 1.20205. . ., CI2 is Clausen's 
function, Cl2(f) = 1.01494..., and 7 = 0.57721... is Euler's constant. The one-loop 
coefficients an and a^^ are similar in magnitude, but the two-loop coefficients bn and by^ 
differ in magnitude by roughly a factor of 7, despite the fact that almost the same number 
of diagrams, with similar structures and magnitudes, contribute. It is also interesting that 
the coefficients in Zjj^ alternate in sign; those in do not. We note that the results above. 



revised relative to our previous analysis [|T2|, are now in complete agreement with the those 
of Ghinculov [|T^, which have also been revised 

Because the decay width F (^H —>■ ff^ is a physical quantity, and all radiative corrections 
that depend only on GpMjj are contained in the factor Zh/Z^ in Eq. (2.6), this factor must 
be finite for e ^ 0. The Z's are finite at one loop, but not at two loops. Hence, the parts 
of the two- loop contributions to Z^ and Z^ that are proportional to 1/e must cancel in the 
ratio. The cancellation is clear if the Z's are written in factored form, 

^ = (1 + aAC + bA'e) (1 + -X'e) + O (A=^) (a = w,H), (2.10) 



and is exact to all orders in A. The complete cancellation of the divergent terms allows us 
to take the limit e — 0. In this limit, — 1 with no pieces left over, and the final ratio is 
independent of the scale n introduced in the process of dimensional regularization. 
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The O {G\Mfj) electroweak corrections to the fermionic decay rates emerge naturally in 
this formalism as the finite ratio 

Zw 1 + auX + huX^ 

This expression for Zh/Z^ automatically resums one-particle-reducible Higgs-boson self- 
energy diagrams in a way that conforms with the standard procedure in Z-boson physics; 
see, e.g., [^. However, it is clear that the resummation contains only limited information 
on higher-order terms. Since we actually have no control of terms beyond O [X^^ , and 
are not aware of a physical principle which would select this as an optimum resummation 
scheme, we expand Eq. and discard terms beyond O (^A^^ = O (Gj^Mfj). This gives 

the alternative representation 

— ^ = 1 + (a^ - aH)X + - hn - a^^au + a^) P (2.12) 
^ 1 + 2.12A- 32.66A2 
^ 1 + 0.013A-0.0013A^ 



The result agrees at O (^Xj with the known one-loop result [[7|,p3| 



- 1 + \ { — -irVs) . (2.13) 



Z^ 8'K^V2 V 2 



III. RESULTS 

A. Limits on perturbation theory 

We are now in a position to explore the phenomenological implications of our results. In 
Fig. 0, we show the leading electroweak corrections to F (^H ff ^ in the one- and two-loop 
approximations with and without resummation of one-particle-reducible higher-order terms 
plotted as functions of Mh- We will concentrate first on the expanded results given in Eq. 
( p.l2| ). While the 0{G']^Mfj) term (upper solid line in Fig. ^ gives a modest increase of 
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the rates, e.g., by 11% at Mh = 1 TeV, the situation changes when the two-loop term is 
included. The importance of this term, which grows as M^, increases with Mh in such a 
way that it cancels the one- loop term completely for Mh = 1114 GeV, and is twice the size 
of the one- loop term, with the opposite sign, for Mh = 1575 GeV. The total correction, 
shown by the lower solid line in Fig. ^ is then negative and has the same magnitude as 
the one- loop correction alone. The perturbation series for the corrections to F {jrl ff^ 
clearly ceases to converge usefully, if at all, for Mh ~ 1100 GeV, or equivalently, for A ~ 10. 
A Higgs boson with a mass larger than about 1100 GeV effectively becomes a strongly inter- 
acting particle. Conversely, Mh must not exceed approximately 1100 GeV if the standard 
electroweak perturbation theory is to be predictive for the decays H //. Note that one 
cannot use the usual unitarization schemes invoked in studies of VF^, Z^, H scattering PJ33 



to restore the predictiveness for the heavy-Higgs width, as no unitarity violation is involved. 

One might expect to improve the perturbative result in the upper range of Mh somewhat 
by resumming the one-particle-reducible contributions to the Higgs-boson wave-function 
renormalization by using Eq. ( |2.1lD rather than Eq. ( [^.121 ). This leads to an increase of the 



one-loop correction (upper dotted line in Fig. while the negative effect of the two- loop 
correction is lessened (lower dotted line) for large values of Mh- However, in the mass 
range below Mh = 1400 GeV, this effect is too small to change our conclusions concerning 
the breakdown of perturbation theory. Moreover, the resummed expression for the one- 
loop terms in the perturbation expansion, when reexpanded to O (G^nM^), does not yield a 
proper estimate for the size of the two-loop terms. There is consequently no reason to favor 
this approach to the present problem. 

It might be argued that the apparent breakdown in the perturbation expansion as judged 
by a comparison of the one- and two-loop terms is an artifact of a small one-loop contribution 
rather than a consequence of large two-loop terms. However, we see no evidence in the 
calculation that there are unusual cancellations in the one- loop corrections. In fact, the one- 
loop contributions an and add in magnitude in the ratio Zh/Z.^}., whereas the two- loop 
contribution is in magnitude subtracted from hn [see Eqs. ( p.9|) and ( p.l2| )]. 



11 



In a previous publication, we looked at other processes to which and Zu contribute, 
e.g., the scattering of Higgs bosons and longitudinally polarized and Z bosons [llO 



In the latter case, even larger coefficients appear in the perturbation expansion when it 
is expressed, as above, in a series in the (running) parameter As = (As/167r^) [^. The 
contributions due to the finite parts of Zu and Z^ are rather insignificant in comparison 
with the finite parts of the unrenormalized two-loop scattering graphs. 

While the factors of (I/IGtt^) occur naturally at each order in perturbation theory, their 
incorporation into the natural parameter A is misleading: the coefficients of the zero-, one-, 
and two-loop terms in the diagonal partial-wave scattering amplitudes, and the one- and two- 
loop terms in the Higgs-boson decay rate calculated above, all have similar magnitudes when 
the series are rewritten as expansions in the physical parameters \ and A = CpM^/ \/2 as 



seen in Eq. ( p.l2|) and the comment [^. It appears, therefore, that A, and not A = A/ (16 



is the natural expansion parameter. 

The high-energy scattering processes give strong evidence for a breakdown of the pertur- 
bation series for a running coupling As(-\/s) ~ 2.3 at either the one-loop or two-loop pO[ 



level. This translates to Mh ~ 380 GeV if the SM is assumed to remain valid for energies ^/s 
up to ~ 5 TeV [0. The Mh upper bound obtained in the present analysis is considerably 
less stringent than the one found in [|l^. However, we emphasize that the result obtained 
here does not depend on the extra assumption about an energy scale; the breakdown of 
perturbation theory is fixed solely by the physical value of Mh- 

B. Comparison with complete one-loop corrections 

The leading corrections discussed so far are independent of the fiavor of the final-state 
fermion. However, from an experimental point of view, they are relevant only for the ti and, 
perhaps, the hh and t^t~ decays of the Higgs boson. It is therefore interesting to compare the 
corrections calculated here with the full one-loop electroweak corrections [^, and the QCD 
corrections that are available for these decay channels P5|,|5B[. In particular, the subleading 
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two-loop electroweak corrections, those of O [G'jpMfjmf) and O (G^mj), are still unknown, 
but one may estimate their likely importance by comparing the top-quark Yukawa-coupling 
correction to the Higgs-coupling correction at one loop. The QCD corrections to the H ^ qq 



modes, where q denotes a quark flavor, are known to 0(a;s) for arbitrary values of nig ||35 



and to O (af) in the limit nig <^ Mh PB|. Their main effect is to replace the pole mass of 



the quark q by its MS mass evaluated at the scale Mh- For completeness, we mention that 
the O (asGpiTii) corrections to the fermionic decay rates are also now available fS^ . 



The results in the preceding section were derived using the Goldstone-boson equivalence 
theorem and neglecting further electroweak corrections of orders g"^, Xg^, etc., as well as 
contributions that involve the Yukawa and QCD couplings of the fermions. Since A, g"^, as, 
and the Yukawa couplings are independent parameters of the theory, our conclusion that 
the perturbation series in A = GfM]j/^/2 fails to converge satisfactorily is independent of 
further corrections involving powers of g and the other independent couplings, though those 
further corrections may be numerically important in applications of the results. 

The use of the equivalence theorem provides a correct framework for calculating the 
leading electroweak corrections — those enhanced by the maximum powers of Mu/My/ — at 
each order. By neglecting subleading corrections that involve g or the Yukawa couplings 
gf = \plnifjv, we expect to obtain a good approximation to the full result provided that 
My/ 1 Mil oc gv/Mii <^ 1 and mj/Mii <^ 1. Because of the high mass of the t quark 
p8| , it is interesting to test the accuracy of the approximation. In Fig. ^, we compare the 
O {GpMfi) correction to F {H ^ ti), already shown in Fig. |^, to the full one- loop electroweak 
correction including the effects of fermions |2^. The full correction was evaluated in the 



on-shell renormalization scheme using rrit = 174 GeV ||3^. We see that the O {GpMfi) 



term underestimates the full one-loop electroweak correction term by 32% (24%) at Mh = 
500 GeV (1 TeV). To check that the difference arises primarily from the inclusion of the 
top quark in the full calculation, and not from the supposedly small contributions from the 
gauge sector of the SM, that is, from a failure of the equivalence theorem, we have carried 
out a complete one-loop calculation using the equivalence theorem with the gauge couplings 
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set to zero p9|, but the top-quark Yukawa coupling retained. The extra contributions are 
O {Gpm'^), and are independent of those considered above. As shown in Fig. the result 
of the calculation reproduces the full one-loop electroweak result very well. The result 
obtained using the equivalence theorem with gt is only 3.9% (1.8%) larger than the full 
electroweak one- loop term at Mh = 500 GeV (1 TeV) for = 174 GeV. The use of the 
equivalence theorem therefore gives a quite accurate approximation to the full theory, even 
for the rather low values of Mh with which we are concerned. The small residual differences 
away from the decay threshold at Mh = ^rrit can be accounted for by the transverse gauge 
couplings, the nonzero masses of the W and Z bosons, and the finite masses and Yukawa 
couplings for the remaining fermions. The extra structure close to the threshold is the result 
of virtual-photon exchange in QED. This generates a Coulomb singularity and a correction 
that behaves near threshold as 1 + aemQt[{'^/'^P) + 0(1)]; where Qt and P are the top-quark 
electric charge and velocity; see left end of the dashed line in Fig. ^. 

In Fig. 1^ we compare the electroweak and QCD corrections. The latter were calculated 
with the asymptotic scale parameter of QCD adjusted to give as{Mz) = 0.118 PD|. The 



one-loop QCD correction to F {H tt) in Fig. ^(a) shows the expected color-Coulomb 
threshold singularity, with a^^Ql replaced in the expression above by (4/3)as(^H)- This 
singularity is associated with the nonrelativistic motion of the quarks. The set of correction 
terms in powers of (27ras/3/5) corresponds to the expansion of a Coulomb wave function at 
zero quark separation. For Mh ^ Sm^, the one-loop QCD correction is negative, with a 
magnitude which increases logarithmically. The two-loop QCD corrections to F {H — > ti) are 
unknown. They are expected to be large close to the tt production threshold, at Mh ^ 2mt. 
At Mh ^ 2mj, the potentially large logarithmic contributions in all higher orders can be 
resummed by using the top-quark MS mass evaluated at the scale Mh-, and the residual 
corrections should be small (see the discussion given below for the hh decay). 

In Fig. ^(b), we repeat the comparison of Fig. ^(a) for the case H ^ hh assuming 
irib = 4.72 GeV |^|. The difference between the full one-loop electroweak correction and the 



O {GpMfj) result is again accounted for at large values of Mh by the omission of top-quark 
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effects in the latter. The spikes in the full correction at Mh = 2Mw and 2Mz originate 
in threshold singularities of the Higgs-boson wave-function renormalization. The dent at 
Mh = 2mt is not accompanied by such a divergence. These features may be understood 
as artifacts of the underlying approximation of treating the unstable Higgs boson as an 
asymptotic state. The QCD corrections are calculated to O (a^) in the MS scheme p5|,|36 



see Eq. (37) of ^11. When the quark pole mass is used as a basic parameter, the largest part 
of the one-loop QCD correction comes from a large logarithmic term — (4«s/vr) In (M///m;,) 
|35| . In general, the large logarithms are of the form (a^/Tr)" ln™'(M///mfc), with n>m. By 
exploiting renormalization-group techniques, these logarithms may be absorbed completely 



into the MS quark mass, mb{fJ,), evaluated at the scale /i = Mh The logarithms are 
resummed to all orders, and the remaining perturbative expansion converges more rapidly. 
The offset seen in the Fig. §(b), with the QCD-corrected decay considerably below the Born 
decay rate, results mainly from the use of mb(M//) instead of mb in the prefactor in Eq. 
( p.7| ). While the effect is large, it is controlled by the resummation. The remaining part of 
the QCD correction at two loops is rather small. 



C. Handling a nonperturbative Higgs 

One would like to be able to describe the Higgs-boson decay to fermions phenomeno- 
logically even in the case of a strongly interacting Higgs sector. In the case of the QCD 
corrections discussed above, the renormalization group provided a physical principle that 
could be used to motivate and organize a resummation of higher-order effects to obtain a 
controlled final expression, even though the corrections could be large when viewed order- 
by-order in as- We do not know of a similar physical organizing principle to use in the 
resummation of higher-order corrections in A, especially as the leading corrections only de- 
pend on one energy scale, namely the mass of the Higgs boson. Any summation of the 
perturbation series will therefore be speculative, and will necessarily be based on the math- 
ematical structure of the series rather a physical argument. We note in this connection that 
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the correction terms in Eq. ( |2.12| ) alternate in sign. This suggests that Pade summation of 
the series might be reasonable. In particular, we can use the two-loop information given in 
Eq. (|2.12| ) to rewrite the perturbative series using a [1,1] Pade approximant that is, as a 
ratio of two first-degree polynomials in A, with the coefficients adjusted to fit the expansion 
in Eq. ( |2.12| ) through order A^. We obtain 

Zw 1 - [{bu, -bn - a-wdh + a'jj)/{a^ - aH)]X 

In Fig. 1^, we compare the Pade-summed correction factor, Eq. (|3.1| ), with the earlier 

results from Fig. ^or Eqs. (|2.11|) and ( p.l2|) . The result suggests that the leading electroweak 

corrections to F (^H ff^ will be quite small even for values of ^1.5 TeV. How far 

the result can be trusted is a matter of speculation. A similar Fade summation of the 

partial- wave scattering amplitudes for W^, Zl, H scattering turns out to give a fairly 

good prediction for the two-loop contribution in terms of the zero- and one-loop terms 

43| , so the method may be more reliable in this rather similar case than our limited input 



information would suggest. If so, the leading electroweak corrections to H ^ ff in powers 
of A or GpMjj will be negligible, when resummed, relative to the corrections introduced 
by the Yukawa couplings and QCD. Only experimental results or reliable nonperturbative 
calculations can resolve this speculation. 

We note in this connection that recent lattice simulations of certain Yukawa models for 
the interaction of the Higgs boson with mirror or reduced staggered fermions suggest that 
the Yukawa couplings cannot be strong, unless the regularization scale is unacceptably low 
PI- 



IV. CONCLUSIONS 

In summary, we have calculated the leading two-loop electroweak corrections to the 
fermionic decay rates of a high- mass Higgs boson in the SM, which are of O {G\Mfj). The 
corrections are negative and exceed the positive O {GpMfj) one- loop corrections in mag- 
nitude for Mh > 1114 GeV. For larger values of M^, the perturbation series is clearly 

16 



unreliable, and the theory becomes effectively strongly interacting. We conclude, given the 
lack of a physical principle that would allow a convincing resummation of the perturbation 
expansion, that a value Mh ~ 1100 GeV has to be considered as a theoretical upper bound 
on Mh beyond which the fermionic decay width of the Higgs cannot be calculated perturba- 
tively. This result is independent of speculations regarding the energy scale up to which the 
SM is valid, as the center-of-mass energy in the Higgs decay is fixed, ^/s = Mh- However, 
there is indication that high-energy interactions in the Higgs sector of the SM can become 
effectively strong, and are not usefully calculable in perturbation theory, for even smaller 
Higgs mass. For scattering processes with ^/s ~ 5 TeV, the critical value for Mh is about 
380 GeV, and for scattering at GUT energies the critical value is less than 160 GeV [1T0[ . 



Clearly, the present-day precision tests of the gauge sector of the SM are not affected by 
such nonperturbative effects. 
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FIGURES 

FIG. 1. The two- loop diagrams that contribute to the wave- function renormahzation constants 



Zh and at O [G'jpMjj) through derivatives of the self-energy functions and n° , Eq. (|]|). 
Heavy (light) lines represent Higgs {w^ and z) bosons. The statistical weights of the diagrams are 



not shown, but may be read off from Eqs. (11) and (12) in |28]. See Appendix A in |^8[ for the 
results. 

FIG. 2. Complete O {OpMl) and O {G^Mfj) correction factors for T {H ^ ff) for 100 GeV 
< Mfj < 1700 GeV. These corrections are universal, i.e., they are independent of the flavor of 
the final-state fermions. In each order, the expanded result given in Eq. ( p. 12 ) is compared to 



the calculation where the one-particle-reducible Higgs-boson self-energy diagrams are resummed 
as shown in Eq. ( |2.11| ). The two- loop correction cancels the one- loop correction at Mh = 1114 GeV 
and is twice as large as the latter, with an opposite sign, at Mjj = 1575 GeV. 

FIG. 3. Comparison of the one- loop results for the ratio F {H tt) /Tb {H tt) obtained in 
various approximations with the full one-loop electroweak result (51, 32, dt^ 56 7^ 0). The solid 
curve (EQT) gives the result obtained using the equivalence theorem with vanishing gauge couplings 
(51 ) 92 7^ 0) and a nonzero top-quark Yukawa coupling corresponding gt to mt = 174 GeV. The 
dot-dashed curve shows the 0(A)=0 {GpMfj) correction from Fig. ||, and it is equivalent to an 
EQT curve with mt=0. 

FIG. 4. Electroweak and QCD correction factors for (a) T (H — > tt), and (b) F [H bb) 
as a function of Mh- universal O [GpMfj) term without resummation (solid line); full one-loop 
electroweak corrections (dashed line); QCD corrections (dot-dashed line); and universal O [GpMjj) 
plus O {G'^pMfj) terms without resummation (solid line). The QCD corrections are evaluated in 
(a) to 0(as) in the on-mass-shell scheme, and in (b) to O (a^) in the MS scheme. The pole-mass 
values mt = 174 GeV and rrif, = 4.72 GeV EO] are used, and the asymptotic scale parameter of 



QCD is adjusted so that Qs(M^) = 0.118 [39|. Note the different scales used in the two plots. 
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FIG. 5. Comparison of the Pade-resummed correction for T[H ^ ff) with the universal 
electroweak correction factors calculated to O (Gi^Adfj) and to O (G^^Af^). 
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